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ABSTRACT 

Nuclear regions of galaxies generally host a mixture of components with different excitation, com- 
position, and kinematics. Derivation of emission line ratios and kinematics could then be misleading, 
if due correction is not made for the limited spatial and spectral resolutions of the observations. The 
aim of this paper is to demonstrate, with application to a long slit spectrum of the Seyfert 2 galaxy 
NGC 1358, how line intensities and velocities, together with modelling and knowledge of the point 
spread function, may be used to resolve the differing structures. In the situation outlined above, 
the observed kinematics differs for different spectral lines. From the observed intensity and velocity 
distributions of a number of spectral lines, and with some reasonable assumptions about structure of 
different subcomponents to diminish the number of free parameters, the true line ratios and velocity 
structures may be deduced. A preliminary solution for the nuclear structure of NGC 1358 is obtained, 
involving a nuclear point source and an emerging outflow of high excitation, ending with shock and 
postshock cloud as revealed by the velocities, as well as a nuclear emission line disk rotating in the 
potential of a stellar bulge and expressing a radial excitation gradient. The method results in a likely 
scenario for the nuclear structure of the Seyfert 2 galaxy NGC 1358. For definitive results an extrap- 
olation of the method to two dimensions combined with the use of integral field spectroscopy will be 
necessary. 

Subject headings: methods: data analysis - galaxies: structure - galaxies: kinematics and dynamics 
- galaxies: active 



1. INTRODUCTION 

The nuclear and circumnuclear activity of galaxies gen- 
erally involves the interplay between a number of differ- 
ent components and phenomena, e.g. a central active 
source surrounded by an absorbing torus, a rotating cen- 
tral bulge, outflowing jets, flowing streams due to the 
action of a bar, or even merging. To separate these dif- 
ferent components and derive their respective line ratios 
and kinematic behaviour is generally difficult due to the 
limited spatial and spectroscopic resolution available. On 
the other hand, such a separation is crucial to the anal- 
ysis of the structures and physical processes involved in 
the nuclear region and their roles in galaxy evolution. 

Then, as to be demonstrated here, such a separation 
could benefit by considering the differences of the veloc- 
ities observed for different spectral lines and be eased by 
models of the activity, smoothed with the point spread 
function (PSF) and fitted to the observations. 

NGC 1358 (Fig. [j) is a barred Sa galaxy hosting an 
active galactic nucleus (AGN) at a heliocentric velocity of 
4028 km s _1 (Theureau et al. 1998), which is included in 
the distance-limited sample of Ulvestad & Wilson (1989) 
as a Seyfert 2. It attracted our attention because of its 
remarkable circumnuclear kinematics which also allows 
demonstrating the utility of the method introduced here. 

Assuming a Hubble constant of 71 km s -1 Mpc -1 
(Spergel et al. 2003), the distance of NGC 1358 amounts 
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to 57 Mpc. Then, 1" corresponds to 275 pc. The outer 
optical disk of NGC 1358, as seen in Fig.[l] including the 
narrow spiral arms, has a radius of 1', or 16.5 kpc. The 
bar along a position angle of 130° has a semi-major axis 
length of 19" (Gerssen et al. 2003) or 5.2 kpc. 

The continuum peak at the nucleus shows a half power 
half width of 1"3. Dumas et al. (2007) derived, by in- 
tegral field spectroscopy, the velocity fields of the ion- 
ized gas in the central regions in the emission line of 
[OIII] A 5007 as well as in the stellar absorption lines 
(the latter displayed only in a poster). The stellar ve- 
locity field shows a kinematic major axis in a position 
angle PA 11.5° with the North side receding. This is 
also close to the assumed PA of the line of nodes of the 
outer disk. The authors derive a symmetric stellar ro- 
tation curve out to ±4" from the nucleus. The [OIII] 
velocity field, on the other hand, shows a kinematic ma- 
jor axis in about PA 90° with the East side receding. 
Narrow band imaging in [OIII] A 5007 shows a slightly 
S-shaped emission line region with a maximum intensity 
1'.'6 from the nucleus, again along a PA of about 130° 
(Dumas et al. 2007). This feature shows clearly deviat- 
ing negative velocities. Dumas et al. (2007) suggest the 
elongated nuclear emission line region to correspond to 
gaseous streaming due to an inner bar inside the Inner 
Lindblad Resonance of the outer bar. 

An HST/STIS spectrum covering the [SII] A 6716, 6731 
lines, within an aperture size of (K'2 centered on the nu- 
cleus, shows double components with the fainter compo- 
nent blueshifted with respect to the stronger (Rice et al. 
2006). 

2. OBSERVATIONS 
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Fig. 1.— NGC 1358 optical image obtained with the EFOSC2 
spectrograph on the ESO 2.2 m telescope. The straight line indi- 
cates the position of the slit (PA 85°) for the spectrum discussed 
in the present paper. The slit crosses the spiral arms at about 35" 
galactocentric distance. 

We have observed NGC 1358 with the EFOSC2 spec- 
trograph on the ESO 2.2 m telescope and obtained a long 
slit spectrum over the nucleus in PA 85° on December 
6th 1996. With this position angle a star was brought 
on the slit, which gave a check on the position angle 
as well as the PSF. The slit cuts the main spiral arms 
about 35" East and West of the nucleus (Fig. [I}. The 
spectrum covered the wavelength region 4075 - 7000 A, 
where the strongest spectral features are the H/3 A 4861, 
[OIII] A 4960, 5007, Ha A 6563, [Nil] A 6548, 6584, 
[SII] A 6716, 6731 emission lines and the Mgb and NaD 
absorption lines (Fig. |2|. In the spatial direction the 
pixel size corresponds to 0'.'26. As flux calibration star 
Hiltncr 600 was chosen. 

Figure [2] shows two fully reduced and flux calibrated 
example spectra obtained from the long slit spectrum. 
To derive the distribution and kinematic parameters, we 
have fitted the galaxy spectrum in two different ways. 
Firstly, by using the standard IRAF routine SPLOT to 
fit Gaussians to the emission and absorption lines, and 
secondly, by using the simultaneous gas and absorption 
line fitting IDL routine developed by Sarzi et al. (2006) 
and Falcon-Barroso et al. (2006). The latter method is 
designed to separate the relative contributions of stellar 
continuum and absorption features and nebular emission 
in spectra of galaxies, while measuring the gas emission 
and kinematics. Accordingly, it uses the routine of Cap- 
pellari & Emsellem (2004) to find the best fitted lin- 
ear combination of single stellar population spectra by 
Vazdekis (1999) to build a stellar spectrum representa- 
tive of the observed galaxy spectrum for this wavelength 
range and with the same spectral resolution. Simulta- 
neously, the nebular emission lines have been fitted, and 
the results yield the gas and stellar kinematics. We found 
that both routines deliver fully comparable gas distribu- 
tion and kinematics, with the advantage that the IDL 
routine has derived a much more robust stellar kinemat- 
ics. 

Figure [| shows the run of [Nil] A 6584, Ha, 
[OIII] A 5007 and H/3 emission as well as NaD absorp- 



tion intensities and velocities over the central region of 
NGC 1358 along the slit as derived with the SPLOT 
routine. The solid line shows the velocities of the stellar 
component derived from all the stellar features. As zero- 
point of the abscissa we have chosen the position of the 
maximum of the continuum (cp. Fig. pj. In Fig. p^the 
emission line velocities, and in particular Ha and]NII], 
show a rapidly increasing rotation curve with a maxi- 
mum about 4" from the center. Dumas et al. (2007) 
show that this nuclear emission line disk should have a 
kinematic major axis directed close to PA 90°, i.e. only 
5° from the PA of the present slit position, but almost 
perpendicular to the kinematic major axis of the stellar 
bulge. 

Within 2" from the center, on the East side, we see a 
reversal, or a dip, in the apparent rotation curve of the 
emission lines. In the NaD absorption line velocities, this 
dip is absent. From Fig. [3] it is clear that the emission 
line velocity dip corresponds to an excess of the emission 
line intensities up to 2" from the nucleus. The NaD lines 
show no such asymmetry and the position of the maxi- 
mum absorption coincides with that of the maximum of 
continuum emission. 

The simplest interpretation of the gas kinematics ob- 
served in Fig. [3] seems to be a central point source and a 
jet-like outflow that are superimposed on a gaseous disk 
rotating in the potential of the stellar bulge. Preliminary 
simulations along the scheme presented below suggested 
that the jet outflow line of sight velocities should be of 
the order of —200 km s . 

A support for this scenario is found in the [SII] emis- 
sion line profiles from HST/STIS spectra published by 
Rice et al. (2006) reproduced in Fig. [4] One spec- 
trum is observed through a l'.'O aperture and another 
through an aperture of 0'.'2. The applied Gaussian fits 
have been constrained to enforce the relative separation 
between the [SII] doublets. In both spectra a strong com- 
ponent is seen in the two lines with very similar flux 
ratios in the two apertures. In the 0'.'2 aperture spec- 
trum a secondary blueshifted component is seen. We 
have assumed 10000 K electron temperature and used the 
atomic parameters of Mendoza (1983) and Osterbrock 
(1989) to calculate the electron density from these spec- 
tra. The calculated relative fluxes yield for the spatially 
extended and dominating component an electron density 
of 793cm -3 . The central component, only appearing in 
the 0?.'2 aperture spectrum, has a density of 1749 cm -3 . 
This denser component is fainter and blueshifted with a 
velocity difference between the strong and faint compo- 
nents of 214 km s -1 , being the innermost part of the jet, 
respectively. 

The stellar velocities show only a very slight variation 
over the region in agreement with the slight positive gra- 
dient expected for the stellar velocity field of Dumas et 
al. (2007) with the present slit inclined 74° to the stellar 
kinematic major axis. The NaD absorption in front of 
the nucleus appears to have a higher redshift than the 
average stellar template redshift, which also is indicated 
by the shape of the line residuals in Fig. [2] 

Most remarkably the velocity dip to the East of the 
nucleus is larger for [OIII] than for the [Nil] and Ha lines 
and apparently still larger for H/3. Outside the dip [OIII] 
shows smaller rotational velocities than Ha and [Nil]. 
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Fig. 2. — Spectrum of the nucleus (top) and 1'.'5 West off-nucleus region of NGC 1358. The best fitted linear combination of single stellar 
populations shows the stellar features fitted simultaneously as Gaussian fits to the emission lines have been applied. The residuals (data - 
total fit) show that the fits have been successful ([OI] A 6300, 6364 have been excluded in this fit). 

rW 



Here we ask ourselves, what additional information on 
the structure and ionization distribution in the nuclear 
region do these differences in velocity behaviour give us? 

3. THEORY 

An explanation for the differences in velocity behaviour 
between the various emission lines would be that the 
whole central region consists of different gaseous com- 
ponents with different velocity behaviour, different exci- 
tation, and different dust absorption, and thus different 
sets of line ratios. When this region is smoothed by the 
PSF, different distortions of the velocity field for the dif- 
ferent lines will result. 

Assume that the region considered covers N pixels and 
contains K separate structural components, each with 

its own line ratios and set of velocities. If and Vj 
are the observed line intensities and velocities at pixel 
number i of a line with wavelength A, we have 
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where JjL is the intrinsic intensity, in the line A, of the 
component k = (1, . . . , K) at the pixel j, and v^j the cor- 
responding velocity, which is supposed independent of A, 
and p(x) the PSF. If L is the number of lines measured, 
we have in the general case KN(L + 1) unknown and 
2NL equations. If the line ratios are all the same in the 
different components, equation pi) is independent of A. 
In the simplest case, where we only consider one compo- 
nent (K — I), we have enough equations to solve for the 
unknown intensities and velocities for each A, which can 
be made by CLEAN-like algorithms (Hogbom 1974). 

In the general case, we do not have a sufficient num- 
ber of equations to solve for the unknown intensities ij^ 
and velocities Vkj. If, however, we increase the num- 
ber of lines measured and at the same time impose some 
conditions on the unknowns in terms of models, we may 
in principle be able to reach a solution. For instance, 
we see from equations ([T]) and ([2| that if the line ratios 

llj 1 ^ /If.j depend on k, i.e., do vary from one compo- 



nent to another, then we have V> XlJ - V^ 2 ' ^ 0. Thus, 
the differences between velocities from different spectral 
lines give us information about differing line ratios be- 
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Fig. 3. — Fluxes and velocities over the central region of NGC 1358 along the slit in PA 85°. The flux scale is arbitrary. The solid curve 
illustrates the velocity of the stellar component as derived from all the stellar features present in the spectra. East is to the right and 
the zero point of the abscissa corresponds to the maximum of the continuum (Fig. [5J , and the NaD absorption line flux values have been 
inverted for presentation purposes. We have derived velocity uncertainties of the order of 10 km s" 1 . 
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4. MODELS 

We will test velocity and intensity observations over 
the nuclear region against a scenario where the dip in 
the rotation curve is caused by an outflow of excited gas 
from a central source at rest with respect to a nuclear 
gaseous disk. We strive to have in this model an absolute 
minimum of free parameters to fit the observations. We 
will first demonstrate the simplest primitive model, here 
called Model 1, to motivate why we introduce our more 
elaborate models. 

4.1. Model 1 



Our Model 1 has the following features: 

a) A stellar bulge with a mass distribution corre- 
sponding to the continuum surface intensity distri- 
bution, giving the gravitational field in the nuclear 
region. As seen in Fig.[5]the observed surface inten- 
sity distribution can be well represented by the sum 
of two exponential disks with scale lengths of 0'.'7 
and 4" respectively and an intensity ratio of 4:1, 
smoothed with the observed PSF. To derive a likely 
potential and rotation curve for this stellar bulge 
we use the relation (2-170) of Binney & Tremaine 
(1987) to give the potential of the sum of two expo- 
nential spheres with scale lengths Qf,'7 and 4'.'0 re- 
spectively, and a mass ratio of 4:1 in conformance 
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Fig. 4.— [SII] A 6716, 6731 spectrum from HST/STIS obser- 
vations of the nucleus of NGC 1358 (Rice et al. 2006), with con- 
strained Gaussian fits. In the bottom panel, the two components 
are separated by 214 km s _1 , and all densities are calculated as- 
suming 10000 K electron temperature. Unconstrained Gaussian 
fits result ip -aljxiQi&t, identical fits. 

with Fig.|5J Although we note that routinely bulges 
are fitted by de Vaucouleurs or Sersic profiles, ex- 
perience has shown that exponential profiles are 
satisfactory to parametrize a bulge (e.g., Fathi & 
Peletier 2003, Balcells et al. 2007), and certainly 
this seems also to be the case for NGC 1358. This 
makes the total mass M x sin 2 («) of the bulge the 
only free parameter of the bulge potential, where i 
is the inclination of the plane of the emission line 
disk to the plane of the sky and we assume that 
the difference in PA between the slit and the line 
of nodes of the disk is small. 

b) In the potential of this bulge there is a rotating 
emission line disk. To avoid edge problems we as- 
sume this disk to be exponential and rotate with 
circular velocities. As the intensities of the disk 
is low compared to the other features, the exact 
scale length is not crucial, and we fix it ad hoc to 
4'/0. Then the only free parameters are the central 
surface line intensities of the emission line disk. 

c) A central emission line point source coinciding with 
the center of the bulge and the emission line disk 
and at rest with respect to that center. 

d) A gas flow, here called the "jet", seen from the 
nucleus out to 7 pixels (IV 8) from the nucleus with 
constant line of sight velocity —214 km s" 1 with 
respect to the point source, which was the velocity 
of the fainter component relative to the stronger as 
seen in the HST/STIS spectrum (Fig. Q). 

The resulting intensities and velocities are convolved 
with the observed PSF as derived from a field star on 
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Fig. 5. — Continuum levels under the emission line Ha fitted 
with the sum of two exponential disks of scale lengths 0"7 and 4" 
and central intensity ratio 4/1, smoothed with the PSF. 
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Fig. 6.— The PSF as derived from the star 150" West of 
NGC 1358. 
the slit (Fig. [6]). 

The line ofsight velocity, or rcdshift, of the central 
source and disk is an open question. Theureau et al. 
(1998) report a redshift of 4028 km s _1 , while Dumas et 
al. (2007) estimate the systemic redshift of the galaxy at 
3997 ± 13 km s~ . Rice et al. (2006) assume a system- 
atic redshift of 4008 km s _1 , and their spectrum (Fig. [5]) 
shows a line of sight velocity of +145 km s _1 for the 
stronger nuclear component with respect to the 4008 
km s _1 systemic velocity, i.e. a nuclear redshift of 4153 
km s _1 for the stronger component and 3939 km s _1 for 
the blueshifted weaker component. We prefer to make 
the systemic velocity of the nuclear disk and the central 
emission source a free parameter in our model. The ad- 
ditional free parameters are the intensities of the central 
source and along the jet for each emission line. 

The fit of this preliminary model is given in Fig. |A1| 
As can be seen, this model cannot reproduce the very 
sharp rise of the emission line velocity curves at the end 
of the jet about 2" to the East of the center. Thus, we 
feel forced to introduce a sharp shock at the end of the 
jet and a postshock cloud. 

4.2. Model 2 
This model has the following features: 

a) A stellar bulge identical with the bulge in Model 1 
presented above, where the total mass M x sin 2 (i) 
of the bulge is the only free parameter. 

b) In the potential of this bulge a rotating exponential 
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emission line disk similar to the one in Model 1, 
again with a scale length of 4'.'0 and the only free 
parameters are the central surface line intensities. 

c) A central emission line point source coinciding with 
the centra of the bulge and the emission line disk 
and at rest with respect to that center. 

d) A gas flow, here called the "jet", seen from the 
nucleus out to a certain number of pixels from the 
nucleus, with constant line of sight velocity —214 
km s _1 with respect to the point source. This again 
is the velocity of the fainter component relative to 
the stronger as seen in the HST/STIS spectrum 
(Fig. EJ). In a chosen pixel at the end of the jet, 
we introduce a sharp shock, with a jump to more 
positive velocities, and in the following pixels, up 
to pixel 8, a postshock cloud with constant velocity. 
The mean velocity in the shock pixel is given some 
intermediate velocity between those of the jet and 
the postshock cloud. 

The resulting intensities and velocities are convolved 
with the observed PSF as derived from a field star on 
the slit (Fig. [6). 

Still, in this way we have tried to keep the free pa- 
rameters to a minimum, and the free parameters of this 
model, to be fitted to the observations, are: 

1. The total mass M x sin 2 (i) of the bulge. 

2. The common redshift of the central point source 
and the center of the emission line disk. 

3. The intensity of the central point source for each 
emission line. 

4. The central intensity of the emission line disk for 
each emission line. 

5. The pixel to contain the shock and the intensity 
and velocity in that pixel. 

6. The intensities along the jet flow for each emission 
line. 

7. The common velocity and the intensities in the 
postshock cloud for each emission line. 

These parameters are varied until a reasonable fit to 
the observed velocities and intensit ies have been ob- 
tained. The fit is displayed in Fig. |A2| The shock is 
placed within pixel 5 with the intermediate velocity of 
— 100 km s^ 1 . This determines the postshock cloud vc- 
locity to about +100 km s^ 1 . As can be seen in Fig. A2 



the line intensities of [OIII] and the velocities in the jet 
are not perfectly fitted, and to place the shock short of 
pixel 5 makes the situation worse. 

Further, it can be seen from the observations that the 
emission lines have long faint wings on the Western (neg- 
ative) side, actually Ha and [Nil] stretching furthest and 
[OIII] having the shortest wing. The run of the Ha and 
[Nil] velocities on the West side show that the main in- 
tensities of this extended emission come from the rotating 
disk. This indicates an asymmetry in the emission line 
disk, and this makes it difficult with the present model to 



fit the velocities away from the center simultaneously on 
both sides, as can be seen in Fig. |A2| Due to the jet, the 
observed velocities on the left side are more influenced 
by the inner part of the disk, and on the right side by 
the outer part. We can then introduce an asymmetry of 
the observations without abandoning the circular sym- 
metry of the disk by describing its surface intensity as 
the sum of two superposed exponential disks with widely 
different scale lengths and varying their central intensi- 
ties independently. This results in differences in the line 
ratios, simulating a gradient of these ratios within the 
disk. This leads us to the next model. 

4.3. Model 3 

This model is the final one which gives us the best fit 
to the observed velocities and intensities in the nuclear 
region of the galaxy. Its features are: 

a) A stellar bulge identical with the bulge in Model 1 
presented above, where the total mass M x sin 2 (i) 
of the bulge is the only free parameter. 

b) In the potential of this bulge there is a rotating 
emission line disk. We assume the intensities in the 
disk to be the sum of two exponential disks with 
differing scale lengths. As the intensities of the 
disk are low compared to the other features, the 
exact scale lengths are not crucial, and we fix them 
ad hoc to 0'.'7 and 4'.'0 (only for historical reasons 
the same as for the bulge). Then the only free 
parameters are the central surface line intensities 
of the two disks (here called Disk 1 and Disk 2) 
representing the inner and outer part respectively 
of the emission line disk. 

c) A central emission line point source coinciding with 
the centra of the bulge and the emission line disk 
and at rest with respect to that center. 

d) A gas flow, here called the "jet" , seen from the nu- 
cleus out to a certain number of pixels from the 
nucleus, with constant line of sight velocity —214 
km s™ 1 with respect to the point source. This again 
is the velocity of the fainter component relative to 
the stronger as seen in the HST/STIS spectrum 
(Fig. [4]). At the end of this jet, at the border be- 
tween two pixels, we introduce a sharp shock, with 
a jump to more positive velocities, and in the fol- 
lowing pixels, up to pixel 7, a postshock cloud with 
constant velocity. 

The resulting intensities and velocities are convolved 
with the observed PSF as derived from the field star on 
the slit (Fig. [6}. 

The number of free parameters of this model is less 
than that of Model 2 and they are: 

1. The total mass M x sin 2 (z) of the bulge. 

2. The common redshift of the central point source 
and the center of the emission line disk. 

3. The intensity of the central point source for each 
emission line. 
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4. The central intensities of the two superimposed 
emission disks for each emission line. 

5. The position of the last pixel of the jet. 

6. The intensities along the jet flow for each emission 
line. 

7. The common velocity and the intensities in the 
postshock cloud for each emission line. 

The fit of Model 3 to the velocities and intensities of 
the various emission lines is displayed in Fig. [7J and the 
parameters are given in Table [I] The table also con- 
tains measurements for the spiral arms with the SPLOT 
routine where they intersect with the present slit. 

5. RESULTS 

Our best fit, Model 3, gives an observed redshift of 
the central source of 4100 km s~\ or 4089 km s" 1 with 
heliocentric velocity correction. The total M x sin 2 (i) 
of the bulge is determined essentially by the velocities 
of Ha and [Nil] at the distance —6" from the center. 
The best fit gives M x sin 2 (i) = 14 x 10 9 M . The 
remaining parameters of this fit are given in Table [l] 
Fig. [8] illustrates Table [T] and gives the intensities and 
velocities of the various features before smearing with 
the PSF. 

The central intensities of the emission line disks are 
essentially given by the run of velocities from —3" to 
-1" and from +2" to +4" from the center. 

The position of the shock falls at the border between 
pixel 5 and 6, i.e. l'/5, from the center. This, in turn, 
forces the postshock velocity up to about +250 km s _1 in 
pixels 6 and 7. 

Figure [7] shows that the velocities and intensities can 
be reproduced with fair accuracy within this nuclear re- 
gion. Finally, the diagnostic line ratios log [Nil] /Ha and 
log [OIII]/Ha, are given in Table [l] for the various com- 
ponents. The H/3 intensities are too weak and uncertain 
to give a unique solution in the present case. 

6. DISCUSSION 

Our "best fit" Model 3 resolves the nuclear structure 
of NGC 1358 into 

1) a central unresolved emission line source (the Nu- 
cleus), 

2) a l'/5 long jet emerging from the nucleus with a line 
of sight velocity of —214 km s~ , 

3) a spherical nuclear stellar bulge containing a rotat- 
ing emission line disk, inclined to the stellar kine- 
matic symmetry plane of the bulge, 

4) an emission line region, the postshock cloud, at the 
end of the jet, and after a drastic jump in velocity. 

With the emission line fluxes for these structures we are 
able to set up the Starburst-AGN diagnostic diagrams 
introduced by Baldwin, Phillips & Terlevich (1981) and 
Veilleux & Osterbrock (1987). Figure [9] shows the diag- 
nostic diagram log [OIII]/H/3 versus log [Nil] /Ha where 
the [OIII]/H/3 ratio is derived from [OIII]/Ha assuming 



large optical depth, referred to as Case B described in 
chapter 4.2 of Osterbrock (1989). Accordingly, for an 
electron temperature of 10000 K, the Ha/H/3 is 2.87, 
which we note is not much different from the case when 
the gas is assumed to be optically thin. As evident, the 
nuclear source and the jet fall in the region of highly 
excited AGNs, while the inner disk and the postshock 
cloud fall close to the region of Low Ionization Nuclear 
Emission line Regions (LINERs). The ratios [Nil] /Ha 
places the outer part of the nuclear disk and the spiral 
arms in the HII region domain. 

The contribution from the underlying Ha and H/3 ab- 
sorption lines only change the derived emission line ra- 
tios marginally. However, the effect of reddening has not 
been taken into account in the emission ratios presented 
in Fig. [9] and Table [T] As we have used the unreddened 
ratio Ha/H/3, our derived [OIII]/H/3 ratios presented in 
the diagram are lower limits in the presence of redden- 
ing. This is an effect that favours our conclusions about 
the nature of the components. 

The nuclear point source shows the highest excitation. 
If H/3 would be included in this analysis, its derived true 
intensity in the Nuclear point source would with Ha give 
the true absorption of the nuclear source of this Seyfert 2 
galaxy. Unfortunately, H/3 is much too faint and uncer- 
tain to be fitted with our procedure. Stronger exposed 
H/3 would be needed. 

The emission line map of Dumas et al. (2007) shows 
that the maximum intensity for our jet lies in about PA 
130°, and there is also a fainter counter-jet in the op- 
posite direction. Thus, our slit avoids the maximum in- 
tensity of the jet and also the counter-jet. As seen from 
Fig. [7J this means that there is no need for a counter- 
jet in our model. However, Disk 1 which is completely 
drowned on he East side by the jet, may in reality corre- 
spond to a trace of an absorbed counter-jet. Its position 
fin Fig. [9 does not contradict this. On the other hand, 
the introduction of a counter-jet in our model would add 
a number of uncertain intensity parameters. 

The HST/STIS spectrum (Fig.]!]) indicates a high den- 
sity for the jet, and the jet intensity maps display a bent 
spiral-like morphology consistent with a precessing jet 
flow. Considering the angle of about 45° between the jet 
axis and our slit, the jet may be surrounded by a more 
widespread outflow, similar to the one observed in the nu- 
cleus of the Seyfert 1.5 galaxy NGC 1365 (Hjelm & Lind- 
blad 1996). Moran et al. (2000) conclude that NGC 1358 
belongs to the class of about 50 % Seyfert 2s that do not 
contain a polarized broad line region (NPBL Sy2). Shu 
et al. (2007) argue that NPBL Sy2s can be interpreted in 
the framework of the AGN unified model if the absorp- 
tion column density is higher at large inclinations and 
the scattering region is absorbed at large inclinations. 
Thus, the inclination of the jet to the line of sight may 
be large and the outflow motion considerably higher than 
214 km s" 1 . 

The emission line disk intensities are completely hid- 
den behind the strong features in the central part and too 
faint to be measured at larger distances from the center. 
They are therefore completely determined by the influ- 
ence of this disk on the observed total velocities. Disks 1 
and 2 represent the inner and outer parts, respectively, 
of this nuclear emission line disk. Their line ratios then 
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Fig. 7. — Best model fit to the emission line intensities and line of sight velocities (Model 3 as described in the text). The dotted line 
indicates the derived rcdshift of the nuclear source. 



TABLE 1 

Parameters for the model presented in Fig. [7] Disk 1 and Disk 2 make up the emission line disk with scale lengths r d of 0"7 

AND 4" RESPECTIVELY (C.F. FlG.[8|. 



Component 


Nucleus 






Jet - 






i — Cloud — > 


Disk 1 


Disk 2 


— Spiral 


arms — > 


Pixel 





1 


2 


3 


4 


5 


6 




7 






128 


-140 


Arcsec 





0.26 


0.52 


0.78 


1.05 


1.31 


1.6 




1.8 


r d = 0.7 


rd = 4 


33 


-37 


Velocity km s — 1 





-214 


-214 


-214 


-214 


-214 


+250 




+250 


Rotating 


Rotating 


-58 


-116 


Flux (Ha) 


4700 


910 


910 


450 


1650 


2100 


1700 




310 


116 


29 


210 


40 


Flux ([Nil]) 


9600 


1400 


1140 


1100 


2700 


5100 


2600 




310 


250 


12 


72 


21 


Flux ([OIII]) 


19400 


3900 


2600 


1720 


3900 


6000 


1500 




600 


160 









log ([Nil] /Ha) 


0.31 






0.28 








0.16 




0.33 


-0.38 


-0.46 


-0.22 


log ([OIII] /Ha) 


0.62 






0.48 








0.02 




0.14 
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Fig. 8. — Illustrating fluxes and velocities of the individual com- 
ponents listed in TablelTJwithout the effect of the PSF. The arrows 
show the fluxes and velocities assigned to the spiral arms which 
are fixed at -37" and 33" galactocentric distance. Filled boxes 
illustrate the nucleus, crosses the jet, and diamonds the cloud, re- 
spectively. 

indicate a rather strong gradient within the nuclear disk 
(see Tabled, or, as mentioned, a difference in line ratios 
between an eventual absorbed counter-jet and Disk 2. 

The disk rotates in the potential of the stellar bulge at 
an angle to the kinematic symmetry plane of the bulge. 
Our slit is oriented close to the kinematic major axis of 
the disk, but close to the kinematic minor axis of the 
stellar bulge (Dumas et al. 2007). The redshift of the 
central source in our model and the systematic redshift 
of the stellar bulge differ with about 100 km s -1 (Fig.[3j). 
The slit crosses the spiral arms close to the photometric 
minor axis of the outer galactic disk (Fig. [lj Their mean 
velocity relative the nucleus is —87 km s _ . In view of 
the fact that our slit is parallel to the line of nodes of 
the emission line velocity field and perpendicular to the 
line of nodes of the stellar velocity field, the difference 
in central velocity between our models and that of the 
stellar bulge, as seen in Fig. [3j may be accounted for if 
our slit misses the central peak passing about 1" to the 
South of the nucleus. This would mean that our derived 
rotation curve is somewhat flatter than the true one and 
the intensities of the central source too faint. The NaD 
absorption lines, measured with SPLOT, show the same 
systematic redshift as the emission lines, which would 



x 



O 




-1.0 -0.5 0.0 
log([NII] ^6548/Ha) 

Fig. 9. — Line ratio diagnostics derived from our ESO spectra. 
The solid curve shows the Starburst-AGN separation lines of Kcw- 
ley et al. (2001), with the shaded region indicating their sample 
galaxies. 

imply that their main absorption comes from the emis- 
sion line disk and not from the stellar bulge. The neg- 
ative mean velocity of the spiral arms speaks somewhat 
against this miscentering view with whatever weight this 
may have. 

7. CONCLUSIONS 

We want to stress that our method, i.a. because of the 
necessity to introduce number of more or less reasonable 
assumptions in order to decrease the number of free pa- 
rameters, does not give a unique solution to the nuclear 
structure of NGC 1358. The differing velocities in the 
different lines give strong constraints for the line ratios. 
Experiments with our model, for example, shows that it 
is the weakness of [OIII] in the postshock cloud and nu- 
clear disk that causes the dip in velocity of the Eastern 
part of the velocity curve to be deeper in [OIII] , and the 
[OIII] velocity curve to be flatter on both sides. But, it 
is necessary to know with good precision the outer wings 
of the PSF in order to account for the influence of the 
strong components on the fainter ones. 

In our particular case the main draw back is that the 
observations are restricted to a single slit. To get all 
necessary information about the velocities and spatial 
distribution of the different components from an isolated 
spectrum is futile, because of the very limited coverage 
and the unknown influence due to the point spread func- 
tion from sources outside the slit. By adding more slits, 
there is still a difficulty to obtain an accurate and com- 
plete coverage (Lindblad et al. 1996). Obviously, the 
ideal is an integral field spectrometer and an extension 
of the method to a two-dimensional treatment, combined 
with comprehensive statistical analysis, covering the en- 
tire region of interest. 

Nevertheless, we argue that, with this analysis, we have 
made a step towards a resolution of the nuclear region of 
NGC 1358 into a number of different components with 
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different velocity behaviour and excitation, where a sum- 
mary of the results are seen in Fig. [9] and Table [T] 
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APPENDIX 

ALTERNATIVE MODELS 

The model presented in this paper have been obtained after a number of iterations fixing the different free parameters. 
Here we show two models based on one single exponential emission line disk, both reproducing the observations to 
different degrees of satisfaction (Fig. Al showing Model 1, and Fig. A2 showing Model 2 as described in the text). 
Although the second model could be seen somewhat comparable to that shown in Fig. [7j we chose Fig. [7] as our 
preferred model since it reproduces most of the observed flux and kinematic features of NGC 1358. 
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Fig. Al. — Fit of Model 1 to the observations. 
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Fig. A2. — Fit of Model 2 to the observations. 



